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Abstract 

We improve on a recent determination of a s from r(T(15) — * j X)/T(T(1S) — * X) with CLEO 
data by taking into account color octet contributions and avoiding any model dependence in the 
extraction. We obtain a s (M x(ls) ) = 0.184^;^|, which corresponds to a s (M z ) = Q.lltf^^. 
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I. INTRODUCTION 



In the early days of QCD heavy quarkonium states (if), bound states of a heavy quark 
and a heavy antiquark, provided an ideal probe of the new theory. Among the interesting 
features, it looked like the strong coupling constant a s could be neatly extracted from the 
ratio T(H — ■> / ygg)/T(H — > ggg), for which both the wave function at the origin and the 



relativistic corrections cancel out 
decays by the Mark II collaboration 



\ The first measurements of J/i/j inclusive radiative 
delivered a photon spectrum not compatible with 
;he early QCD predictions. To a lesser extent, this was also the case for bottomonium (see 
5| and references therein). With the advent of Non- Relativistic QCD (NRQCD) 0,3, ^ 
was understood that color octet contributions, which were ignored in the early approaches, 
become very important in the upper end-point region of the spectrum 8|. When the color 
octet contributions are properly taken into account, a very good description of the photon 
spectrum can be obtained from QCD, at least for the Y(IS') state [9j]. Color octet contri- 
butions also affect the ratio T(H — > / ~fgg)/T(H —>■ ggg) and are parametrically of the same 



order of the re 
to be small 



ativistic corrections. They have so far either been ignored [h]] or estimated 



llj in the available extractions of a s from this ratio. In this paper we take into 



account recent determinations of the T(IS') color octet matrix elements both on the lattice 



12H and in the continuum 



13] , which indicate that their contribution is actually not small. 



This, together with the good theoretical description of the photon spectrum [9|, allows for 
a consistent extraction of ck s (Mx(is)) at NLO in the NRQCD velocity counting. We obtain 
a precise determination of it from the recent CLEO data [lo| . 



II. a s EXTRACTION 

The experimental value of i? 7 = r(T(lS') — > 7A)/r(T(lS') — > X), X being hadrons, 
has been determined most recently in [lO] 2 . We will use only the value obtained from the 
Garcia-Soto (GS) parameterization of data [9], which follows from a QCD calculation. This 

1 See (sj for an update on a s extractions from heavy quarkonium systems. 

2 7 stands for direct photons only and the contributions T(IS') — > 7* — > X have been subtracted in the 
denominator. 
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IS 



R e * p 



0.0245 ±0.0001 ±0.0013, 



(1) 



where the first error is statistic and the second systematic. 
Our starting point is the expression: 

r(T(15)^ 7 X) 36 e 2 b aN 



1?- 



T(T(1S)^X) ~ 5 a s D' 



N = l + C ggi - + C 



71 



•'Pi( 3 Si)^--Pi(3Si) 



71 



+ —a 



7 8 ( 1 So)^0 8 ( 1 5o) 



7T 



H — C 7 o 8 (3p )7?.08( 3 Po) + CjvC'u ), 

C 999~ + ^i( 3 5i)^-Pi( 3 5i) + — C08( 3 5i)^0 8 ( 3 5i) 
7T Q^g 

7T 

±— Co 8 ( 3 Po)^0 8 ( 3 P ) ± Cd(w ), 



+ -q 

a, 



(2) 



(3) 



(4) 



where n/ = 4 is the number of active flavors, a the fine structure constant, e& = —1/3 the 
bottom quark electromagnetic charge, a s = a s (Mx(is)) is the strong coupling constant cal- 
culated at the T{1S) mass, M T (is) = 9.46 GeV, C Pl( 3 Sl ) = -(19vr 2 - 132l/(12vr 2 - 108) 

T|, C7O 8 ( 3 p ) = 189/(4tt 2 - 36) Q, Co 8( 3 Sl) = 
fcl = 81/(8tt 2 - 72) (l6|, C 08(3Po) = 567/(8vr 2 - 72) U, 
13, U, C 5!?!? = 3.79 ± 0.54 (for n f = 4) [7J, U, ft = 



CWft) = 27/(47r 2 - 36) 
81n / /(207r 2 - 180) 16], C 0g( i 



.'/.'/I 



—1.71 (for nj 



4) 



(T(lS)|0|T(lS))/(mf d (T(15)|Oi( 3 Si)|T(lS))), where Ad is the difference in dimension 
between the operators O and OiQSi), m& is the bottom mass, and the (T(1S')|0|T(1S')) are 
NRQCD decay matrix elements [7\. If we adopt the counting of {?| and a s /7T ~ v 2 , then the 
expansions ([3]) and (HI) are valid up to order v 2 . On(v 3 ) and Cd(w 3 ) account for higher-order 
corrections of order v 3 . In the following we will assume v 2 = 0.08. 

In order to obtain a sensible extraction, the ratios of NRQCD matrix elements 1Z must 



be correctly estimated. 7^(3^) can be related to the binding energy [l9j, |20J. 1Z 



K 



O s ( 3 Si) 



have been ca 



and 



culated on the lattice [12|. T^o B ^s ) an d 7^o 8 ( 3 Po) have been estimated 
in the continuum [13j. The continuum calculation and one of the lattice calculations of 
TZo^So) are compatible. We will present two different extractions: C (for continuum) and L 
(for lattice). Extraction C uses all the weak-coupling expressions available, in the same way 
they were used for the description of the photon spectrum [9J, and the lattice calculation 
of 12J for 7?-o 8 ( 3 5i)- Extraction L uses all the lattice calculations available, and NRQCD 
velocity scaling to estimate TZo s (sp y In both extractions, we do not expand the 0{y 2 ) 
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terms in D: it turns out that even though they are individually small they add up to give 
a contribution comparable to one; we will comment on this in Sec. IIHI 



A. Extraction L (for lattice) 

Concerning the ratios TZo, we will take them in the following ranges: 

< n 08{ i So) < 4.8 x 10- 3 , (5) 

< ^ 08( 3 5l) < 1.6 x 10- 4 , (6) 

-2.4 x KT 4 < ^o 8( 3 Po) < 2.4 x 10~ 4 , (7) 

-0.16 < K Pl[ 3 Sl) < 0. (8) 

quations ([5]) and ([6]) correspond to the maximum values obtained in the lattice calculations 



12J taken with a 100% uncertainty. Equation ([7]) follows from the naive counting 



|^O 8 (3P )l - g 



(r /b |0 8 ( 1 P 1 )l^) 



i 



1 V 

- x , (9) 

9 2N r V ; 



m^T(15)|0 1 (35 1 )|T(lS)) 
taken with a 100% uncertainty. The first equality is due to spin symmetry and is valid at 
leading order in the velocity expansion, in the second one we have evidenced the color factor 
1/(2JV C ) [12] (iV c = 3 is the number of colors). 

Equation (jS]) follows from the Gremm-Kapustin relation [19I . I20I ] in the weak-coupling 
regime, 

n VieSl) = ^ « -v\ (io) 

taken with a 100% uncertainty, i^bin stands for the binding energy. The operators On{v 3 ) 
and Od{v 3 ) are taken in the range 

- 0.04 < O n {v 3 ),O d {v 3 ) < 0.04, (11) 

which encompasses both 0(v 3 ) and 0(af) corrections. 
The fine structure constant is taken at the T(1S) mass 

«(^T(is)) = (12) 

We evaluate the theoretical side of Eq. (j2j) as it stands, without further expansions. 
Taking the central values for C ggg and in Eqs. ©-(JED, dlT]) and (pT|) , we obtain 

a s (M T(ls) ) = 0.1885 . (13) 
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The uncertainty on a s induced by a given parameter is evaluated by varying it in the range 
and keeping all other parameters at their central values. We obtain 



Os( l s ) 


= 0.0047, 


'7? i tt s 


= 0.0019, 


^0 8 (3P 0) «S 


= 0.0032, 




= 0.0106, 


1-77 , 


= 0.0041, 


%(» 3 ) s 


= 0.0031, 



6 Cggg a s = 0.0025, (14) 

(15) 
(16) 
(17) 
(18) 
(19) 
(20) 

5 R c r a s = 0.0089. (21) 
We sum up linearly the errors 8n . and 5n , , which are correlated, and then all the 

°8( S o) °8( S l) 

errors quadratically, obtaining 

a s (M r{ls) ) =0.189 ±0.017. (22) 

The dominant error comes from the uncertainty in lZ>p 1 (3s 1 )- We can reduce this uncer- 
tainty, by noticing that for H Vl (3 Sl ) we have an explicit expression, Eq. ffTOl) . that we have 
only partially exploited. Indeed, in the weak-coupling regime, the exact form of is 
known. At the order we are interested in, it holds that 

£bin (C F a s ) 2 4 

: > °f - -, l/oj 

m b 4 3 

where a s is evaluated at the scale Mrfis)Cpa B /2, the typical momentum-transfer scale in a 
Coulombic bound state. From Eq. ( 1221) . we obtain: 



a s (M T(ls) C F a s /2) = 0.311 ± 0.032 , (24) 

which gives 

n Pl( s Sl) = ^ = -(0.0431°;™). (25) 
nib 

Using this value for TZ-p 1 ^s 1 ) an d performing again the above calculation we obtain the new 
central value 

« S (M T(15) ) = 0.183, (26) 
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and the new uncertainties 

6 Cggg a s = 0.0026, (27) 

(28) 
(29) 
(30) 
(31) 
(32) 
(33) 

^a s = 0.0085. (34) 
Summing up the errors like before, we obtain as our best estimate 

a s (M T(ls) ) = 0.183 ± 0.013 . (35) 
This corresponds to a strong coupling constant at the Mz mass of 

a s (M z ) = 0.119 ±0.005. (36) 

B. Extraction C (for continuum) 



K„ , a s 
O 8 (!s ) 


= 0.0040, 


^0 8 (3 Sl ) a s 


= 0.0026, 


^O 8 (3p ) a s 


= 0.0027, 


1? Ola 

St 3 Si) S 


= 0.0014, 




= 0.0044, 


%(» 3 ) S 


= 0.0033, 



In a weak coupling analysis, (Y(lS')|0 1 ( 3 5'i)|Y(lS')) can be calcu 
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ated in perturbation 
3 . In order to follow 



theory of a s (rribv). A NNLO expression is necessary at 0(v 2 

the same procedure as in j^] , we multiply the leading order term in the decay widths by the 
NNLO expression for (T(1S')|Oi( 3 5'i)|T(1S')) and the a s correction to the decay widths by 
the LO expression for that matrix element. If we factor out the NNLO matrix element, this 
produces a shift iV -> N + 5N and D -> D + SD in <^ and flU): 

5N = C ggj ^5, (37) 

7T 

5D = C 999 ^5, (38) 

7T 



with 

. (T(lS)\0 1 ( 3 S 1 )\T(lS)) LO 



(T(15 , )|0 1 (35 , 1 )|T(15')) NNLO 



1 . (39) 



We count a s at the soft scale as order v. 
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For the central values of the objects below we take exactly the same ones used in [9[, namely 

5 = -0.57, (40) 

n Vl{3Sl) = -0.015, (41) 

n o ^s ) = 0.0012, (42) 

n Os( t P0) = 0.0011. (43) 



For TZo^s!), w e use the hybrid algorithm output of the lattice calculation [12j] 4 . 

^ 08( 3 5l) = 8 x 10- 5 . (44) 

Using those values, we obtain 

a s (M T{15) ) = 0.185. (45) 

In order to associate errors to these central values, we move the values of the objects 
below in the following ranges: 

0.18 < a s {m b v) < 0.38, (46) 
0.32 < a s (m b v 2 ) < 1.3, (47) 
0< ft 08(35l) <1.6xl0~ 4 . (48) 

The wide variation range of a s (rribv) and a s {m b v 2 ) is expected to account for O(Aqcd) 
uncertainties in the weak coupling estimates of O 8 ( 1 5 , ) and O 8 ( 3 P )- The upper limit of 



Os( 3 Si) corresponds to twice the largest value obtained using the lattice algorithms in [12 ] . 

The uncertainty on a s induced by a given parameter is evaluated by varying it in the 
range and keeping all other parameters at their central values. We obtain 

8 Cggg a s = 0.0009, (49) 

A rv —+0.0006 (cn\ 

y a a (m b u)«s — -0.0064' \ ov J 

A o _+0.0083 /ri\ 

<Ja s (m b v 2 ) u; s — -0.0076; \ OL ) 

5 n , a s = 0.0016, (52) 

X n _+0.0035 (ro\ 

° n o N (v3) as ~ -0.0034 > l M J 

X r, _+0.0026 (r A \ 

° n o D { V 3) as ~ -0.0025 > \P*) 

<Wa s = 0.01. (55) 



The hybrid algorithm is selected because compares well with the continuous estimate for 7?.o s ( 1 s )- 
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We assume these errors to be independent and sum them up quadratically, obtaining 

« s (M T(ls) )=0.185i^. (56) 
This corresponds to a strong coupling constant at the M z mass of 

a s (M z ) = 0.1201°;^ . (57) 

III. DISCUSSION 

We have presented two extractions of a s at NLO in the NRQCD velocity counting, the 
main differences being the values assigned to the NRQCD matrix elements. The two out- 
comes are very close so we take the average as our final central value. Since the two extrac- 
tions are not completely independent we take as our error the range of the two determina- 
tions. Then our final value is 

a a (M T(15) ) =0.1841°;°^, (58) 

which corresponds to 

a s (M z ) = 0.1191°;™?, (59) 



which is very close to the central value of the PDG [23| with competitive errors. The 
key ingredients to get these numbers are the precise CLEO data [lOj], the use of a QCD 
calculation (called GS model in 1Q| ) to extrapolate the photon spectrum at low z, and 
accurate estimates of the color octet matrix elements, which have been possible thanks to 
recent lattice and continuum estimates. Concerning the matrix elements, our results are 
rather insensitive to the values of 8 ( 1 S'o) and 8 ( 3 -Po) (note that the upper limit given in 
Eq. (j7j) for TZo 8 (3p ), based on the scaling (jUj), is smaller by a factor five than the continuum 
estimate (j4"5|) ). but would be sensitive to large values of Os( 3 Si). However, the lattice values 
for Os( 3 5'i), which we have used, turn out to be much smaller than what NRQCD velocity 
scaling rules suggest, and do not have a major impact in our results. 

How reliable is our extraction? Our determination is valid at next-to-leading order in 
a s (mb) and in v 2 . At this order, terms corresponding to new qualitative features appear 
(radiative, relativistic, octet corrections), each of them of natural size, but whose sum is 
of order one and hence large. This is not unusual. It is crucial, however, that higher- 
order corrections, those that we have generically labeled as On and Od, are small. This 
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is expected because higher-order corrections do not introduce new qualitative features. In 
Ref. 24J], higher-order corrections in the velocity expansion can be found for Od- These 
are not the complete set of corrections entering in 0£>, since higher-order a s corrections to 
the lowest operators are missing. Anyway, using the analogous of Eq. ffTU]) to estimate the 
matrix elements and taking E^ in /mb « —0.04, the corrections calculated in [24j amount to 
about 0.02, which is consistent with the estimate (fill) . Analogous corrections for On are not 
known. At present, the main uncertainty in our extraction of a s comes from the systematic 
uncertainties in R: xp - 



Let us next compare our extraction to two previous related ones [lOl. 111]. 
Concerning the extraction of the CLEO paper 



lOj^here are two main differences, (i) 
18j], in which the NRQCD color octet 



On the theoretical side an old formula was used there 
operators were ignored. This introduces large theoretical uncertainties. In practice, however, 
we find that numerically they are not so important for the final result, (ii) For the total 
radiative width, two numbers are quoted depending on whether the so called Field model 
25| or the GS model, which is in fact a QCD calculation, are used for the extrapolation 
of the photon spectrum at low z. The final number is given as the average of the two 
procedures. We believe that the use of the Field model, which uses a parton shower Monte 
Carlo technique to incorporate the effects of gluon radiation by the outgoing gluons in the 
decay, introduces an unnecessary model dependence that moves the actual central value and 
artificially increases the errors. Our final results are similar to the ones presented in 10j for 
the GS model. 



Concerning the extraction of 111] , which is used by the PDG 23j], there are three main 
differences, (i) On the theoretical side, the color octet NRQCD matrix elements are ignored 
in r(T(15') — > jX), whereas we find that they contribute between 30% and 100%. 5 (ii) 
Older data are used, which are fully consistent with, but not as precise as, the more recent 
ones, and an older analysis 261 ]. which relies on the Field model for extrapolations to low z. 



(iii) The extraction is actual 



y done from T(T(15) -»■ 7 A)/r(T(15) -> We believe, 



contrary to a statement in ll| , that the latter increases rather than decreases the theoretical 



uncertainties associated to color octet operators. Indeed, whereas the ratio radiative/total 



5 The color octet contributions in r(T(15) — > X) are estimated to be 9%, whereas ours turns out to be 
between 50% and 160%. 
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has the same color octet operators in the numerator and denominator except for one, the 
ratio radiative/leptonic (total/leptonic) has two (three) different color octet operators in the 
numerator and denominator. Furthermore, the leptonic width is known to suffer from large 
higher order corrections in a s (see 27J] for a recent discussion), which introduces further 
uncertainties. 



IV. CONCLUSIONS 

We have improved on current determinations of a s (M r ^ ls ^) from radiative decays of T(1S) 
by avoiding any model dependence and by taking into account recent estimates of color octet 
operators. The value we obtain, a s (Mz) = 0.119^o;qq5, is close to the PDG average with 
competitive errors. 
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